Initial static susceptibilities 
of nonuniform and random Ising chains 

Oleg Derzhko, Oles' Zaburannyi 

Institute for Condensed Matter Physics, 
1 Svientsitskii St., L'viv-11, 290011, Ukraine, 

J. W. Tucker 

Department of Physics, The University of Sheffield, 
Sheffield, S3 IRE, United Kingdom 

February 7, 2008 

Abstract 

Within the conventional framework of standard linear response theory we have derived 
exact results for the initial static susceptibilities of nonuniform spin-^ Ising chains. The 
results obtained permit one to study regularly alternating-bond and random-bond Ising 
chains. The influence of several types of nonuniformity and disorder on the temperature 
dependence of the initial longitudinal and transverse static susceptibilities is discussed. 

PACS codes: 75.10.-b 

Keywords: Regular alternating-bond Ising chain; Random-bond Ising chain; Susceptibility 
Postal addresses: 

Dr. Oleg Derzhko (corresponding author) 
Oles' Zaburannyi 

Institute for Condensed Matter Physics 

1 Svientsitskii St., L'viv-11, 290011, Ukraine 

Tel: (0322) 42 74 39 

Fax: (0322) 76 19 78 

E-mail: derzhko@icmp.lviv.ua 



Dr. J. W. Tucker 

Department of Physics, The University of Sheffield 

Sheffield, S3 7RH, United Kingdom 

Tel: +44 114 222 4273 

Fax: +44 114 272 8079 

E-mail: j. w. tucker© Sheffield, ac. uk 



1 



1 Introduction 

During the last 40 years considerable effort has been devoted to the study of random mag- 
netic systems. In order to obtain a theoretical understanding of various phenomena that can 
occur due to the presence of disorder it is instructive to examine simple models that allow one 
to calculate observable quantities exactly. Among such models, perhaps the simplest is the 
nonuniform spin-^ Ising chain specified by the Hamiltonian 

n 

where the site index n of the summation takes values 1 to AT — 1, there being a total of N lattice 
sites. J n is the exchange coupling between the spins at sites n and n + 1, and s x has eigenval- 
ues ±1/2. The partition function of this model is readily calculated, and for open boundary 
conditions is [0J Z = Tr [exp(—H /kT)] = 2 IIn[2 cosh( J n /4kT)], which immediately yields 
the entropy S = k [In 2 + £ n In 2 cosh(J ft /4fcT) - £n(J n /4fcT) tanh( J n /4kT)), the specific heat 
C = kJ2 n [(Jn/^kT)/ cosh( J n /4kT)} 2 , and the correlation functions 

Here, the angular brackets ((. . .))o denote (l/Z )Tr [exp(H /kT)(. . .)]. Regarding the exchange 
couplings in Eq. (1) as random variables having a probability distribution p(Ji, . . . , Jjv-i)> 
one is able to study the random-averaged quantities S, C, or (s^Sp)o, defined through Q = 
/ dJi . . . dJjv_! p( Ji, . . . , Jat_i)Q. In the following, attention will focus on a study of the initial 
(zero- field) static susceptibilities of nonuniform and random Ising chains described by Eq. (1). 
In the presence of a small external field h a , a = x,y, z, static susceptibilities may be defined by 
Xf3a(h a ) = dmp/dha, mp = (1/N) J2 n ( s n) (i n um ts where the product of the g-factor and the 
Bohr magneton is unity). The angular brackets ((...)) denote thermal averages with respect to 
the total Hamiltonian H = Ho + H a , H a being the Zeeman term —h a J2n s n- The initial static 
susceptibilities are the limits X/3a(h a )]h a -*0j to be denoted by X/3a- 
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The static susceptibilities of random Ising chains have been studied by a number of authors. 
The Ising chain with random exchange coupling in a longitudinal field was considered by Fan 
and McCoy 0. Assuming that each exchange coupling was an independent random variable 
having a probability density function with narrow width, Fan and McCoy were able to study 
the influence of the randomness on the longitudinal static susceptibility. The ground-state 
transverse static susceptibility for the random-bond Ising chain was examined by Barouch and 
McCoy [[| . Zaitsev studied the transverse static susceptibility for the random-bond Ising chain 
for various types of disorder [[|. The thermodynamics and spin correlations of a model of 
solid solution based on an Ising chain were considered in Ref. ||. Some results for the initial 
longitudinal and transverse static susceptibilities of the random Ising chain in a transverse field 
were collected in Ref. . Thermodynamic and dynamic properties of a random Ising chain in a 
transverse field for arbitrary disorder may also be examined numerically 0. Recently, Idogaki, 
Rikitoku and Tucker M have calculated exactly the initial transverse static susceptibility for 
the nonuniform model of Eq. (1), and for the random version of this model in which two types 
of exchange couplings J\ and J2 are randomly distributed throughout the chain. 

The object of our paper is two- fold. First, in Ref. || the derivation of the initial static 
transverse susceptibility was achieved by a somewhat unconventional technique based on an 
accurate extraction of the linear field- dependent term from a Callen-Suzuki identity for the 
magnetization. Here, we show how the result may be obtained more naturally within the 
framework of conventional linear response theory in which the susceptibility is related directly 
to the double-time correlation functions. In addition all components of the susceptibility tensor 
will be considered (Section 2). The second objective of this work is to apply the theory to 
chains having various types of disorder, to identify the extent to which qualitative features of 
the susceptibility are influenced by the character of the disorder, and to compare the influence 
of regular alternation and random disorder. Unlike, the bimodial distributions of disorder, 

3 



hitherto studied ||, more general distributions of disorder cannot be treated analytically, but 
require numerical computation. The result of such computations for chains having Gaussian 
and Lorentzian disorder are reported for the first time, in Section 3. 



2 Initial static susceptibilities 



of nonuniform Ising chains 



Our derivation of the initial static susceptibilities of a nonuniform Ising chain employs as its 



starting point the well-known expression 



XPa{h a )- kTN Y. 



4{t) ee e^e"^, 



(3) 



relating static susceptibilities to spin correlation functions. For the initial static susceptibilities 
(h a —>■ 0) it follows, from Eq. (3), that 



Xxx — ^rp 5Z(' S n' S p)0) 



(4) 



and 



Xz 



1 1 



S/o ldr( ^(-^)^ 



JtH Q J3 -itH 



(5) 



since (s") = 0. 

The integral in Eq. (5) may be effected by employing the Van der Waerden identity 
exp(Xs x ) = cosh(A/2) + 2 sinh(A/2)s :r . Using this identity, one finds that 



= COS [t (Jj-isj.! + JjS* +1 ) 



Sj + sin 



t ^J 7 _isJ_ 1 + JjSj +1 



J 3 ' 



(6) 
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which, since cos(Xs x ) = cos(A/2) and sin(As a ) = 2 sin(X/2)s x , gives 

m = 



cos 









- cos — - 


2 


2 




-l £<A' 




— cos — - 


2 


2 



- A sin — sm -^s J _ 1 sJ +1 j s~ 



+2 ( sin ""^ cos -rrs x _ 



tJj—i . tJj , „ 



(7) 



With the aid of this expression, the integration in Eq. (5) is readily perform to yield the desired 
result 

J-n, — 1 H~ J<n 



Xz 



( sinh n 'l T Jn sinh Jn ~ 1+ Jn 



+ 



2kT 



n,p 



\ Jn—1 Jn Jn—1 J-n 



\ S n S p)o 



\ Jn—1 Jn Jn—1 Ji 



■ Jn-l+Jn \ 

2kT \ / Q x z x z\ 



-2i 



/ 1^ Jn-l-Jn _ 1 |^ Jn-l + Jn _ 1 \ 

H t ! T \ S n-l S n S p/0 



\ Jn—1 Jn 



+2i 



/ cosh 



Jri—1 J-n 

2kT 



1 cosh 



Jn—1 ~t~ '-'n 
2fcT 



1 



( S n S n+l S p)o 



\ Jn—1 Jn Jn—1 Jn 

Since H does not contain s y or s z , and s y s z = (i/2)s x , Trs a = and {s a ) 2 = 1/4, it is readily 
seen on using a representation in which s x is diagonal that the correlation functions appearing in 
Eq. (8) can be simplified. Namely, {s z n s z p ) = (1/A)5 np , «_i« + iSpo = (1/A)5 np (s x _ 1 s x +1 ) , 
( s n-i s n s p)o = (i/ 2 )M<_i<)o, ( s n s n+i s p)o = (i/2)5 np (sX+i)o- The pair correlation func- 
tions in s x that remain may then be evaluted using Eq. (2). In this way, it follows after 
straightforward calculations, that the initial longitudinal and transverse static susceptibilities, 
Eq. (4) and Eq. (8), of the nonuniform Ising chain may be expressed in the form 



AkT 



n p>n 



tanh " ^ ( — tanh . . . ( — tanh 1 

AkT J \ AkT J V AkT 



Xz 



1 y, ( tanh ^ ~ tanh ^ | tanh ^ + tanh ^ 



(9) 
(10) 



2iV n y J n _i — J n Jn-l + Jn J 

The result, Eq. (10), for Xzz coincides with the result obtained by Idogaki, Rikitoku and Tucker 
H using a different approach. 

The off-diagonal components Xa/3 can be shown to be zero as follows, xpa contains 
Tr[exp(— Ho/k^s^^s }. Again, using a representation in which s x is diagonal, this trace 
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is zero if /3 — x and a = y or z. Thus \ xy and \ xz are zero. Also, from the definition of the 
interaction representation, s(t), and the fact that a trace is invarient under cyclic rotation of 
the operators under it, it follows that Xyx and Xzx are also zero. Further, Xzy is given by Eq. 
(8) with Sp replaced by s y p . All the correlation functions on the right hand side then reduce to 
zero, so Xzy = 0. Finally, Xyz = Xzy (= 0) and Xyy — Xzz from the symmetry transformation 



,1/ 



-s x , s yl = s z , s zl = s y . 



One immediately observes for the regular chain, J\ = J% = . . . , that Eq. (9) gives 

1 



X Jl 

A.xx 



AkT 



2 / 7i x q 

l + -Y(N -q) f-tanh4^ 



AkT 



e 2kT 



1 + tanh / 4A;T 



\_ 2 tanh ^ 



in the thermodynamic limit. For the regular alternating-bond chain, J% = J3 
. . . , Eq. (9) yields 



Ji — Ja 



v JlJ2 

An 



1 (l -tanh ^l.) (l -tanh ^ 



4A;T 



1 - tanh tanh ^ 



(12) 



and for the case of the regular alternating-bond chain J x = J 2 = J5 = J§ . . . , J3 = J4 = J-j 



j\ Ji J3 J3 



4A;T 



1 - 



tanh g|g + tanh ^ 
1 - tanh ^ tanh 



I (tanh ^ + tanh ^) +2 tanh 2 ^ tanh 2 ^ 



1 - tanh 2 ^ tanh 2 ^ 



(13) 



Further, for the transverse components, Eq. (10) yields the results: xll H3> xil J2 [PL xil JlJ:iJ3 



4Az2: ~r 4A22 ~ 2^-^^ ' Azz 3A22 ~r 3A2Z ' 3A.ZZ clUl 

It is worthwhile to note, that Xzz, Eq. (10), does not depend on the sign of the intersite 
coupling, whereas Xxx, Eq. (9), does. Thus, according to Eq. (11) for the regular chain, one 
finds immediately that for ferromagnetic coupling (Ji < 0) Xxx — > °o as T — > 0, whereas for 
antiferromagnetic coupling ( J x > 0) Xxx -> as T -> 0. In the case of the regular alternating- 
bond chain J\ — J 3 = . . . , J 2 — Ja — . . . , Eq. (12) implies that Xxx diverges as T — ► 
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if both Ji and J 2 are ferromagnetic, whereas if either, or both, of the intersite couplings is 
antiferromagnetic, it does not diverge. From Eq. (13) for the alternating-bond chain J\ = 
■h = -h — Jq — ■ ■ ■ , -h — Ja = -h — Js — ■ ■ ■ , one finds that if both, or either, of J\ 
and J3 are ferromagnetic Xxx diverges as T — > 0, whereas, it does not, if both couplings are 
antiferromagnetic. These results are to be expected if one considers the ground state of these 
regular alternating-bond spin chains. For example, it is clear that the ground state of the chain 
Ji = J3 = . . . , J2 = J4 = ■ ■ ■ , with Ji < 0, J2 > 0, is of an antiferromagnetic type, whereas 
the ground state of the chain J x = J 2 = J5 = Jq = . . . , J3 = J4 = J7 = Js = ■ ■ ■ , having 
Ji < 0, J 3 > 0, is ferromagnetic in character. 

3 Initial static susceptibilities 
of random Ising chains 

Consider an Ising chain in which the exchange couplings are random variables evenly distributed 
with a probability p(Ji, . . . , J/v-i) = Y\ n p{Jn)- The random-averaged susceptibilities are given 



by 




(14) 



(in the first equality use has been made of the fact that — tanh( J p /4kT) is the same for all sites 



p), and 




) 



+ 



( 



tanh + tanh -f^ 



Ji + J 2 



) 



(15) 



Besides the probability distribution 



p{J n ) = c5(J n - Ji) + (1 - c)5(J n - J 2 ), < c < 1 



(16) 



7 



for which Xzz was examined in Ref. ||, the Gaussian distribution 



P(Jn) 



1 



exp 



(Jn ~ JO? 

2a 2 



(17) 



and the Lorentzian distribution 



P(Jn) 



i r 



(18) 



vr (J n - Jo) 2 + T 2 



centred at Jo with a strength of disorder controlled by a 2 and T, respectively, will be considered. 

Consider first the longitudinal susceptibility Xxx, Eq. (14), (Figs. 1-3). The low-temperature 
behaviour, discussed in Section 2, for the uniform chain arose because tanh(J/4/cT) — > sgnJ as 
T — > 0. As a result, in this limit Xxx diverges more rapidly than 1/T for J < 0, but does not 
diverge for J > 0. For the random case, the essential quantity influencing the initial longitudinal 
static susceptibility is | tanh(J/4A;T) |. If there are present, exchange interactions of opposite 
sign, then this quantity is < 1 as T — ► and a qualitative change in the temperature dependence 
of Xxx from that of the uniform chain results. The low temperature behaviour is now of the 
paramagnetic type, 1/T, (the dotted lines in Figs, lb, lc and 2). A similar qualitative change 
in the temperature dependence occurs in the presence of randomness defined by Eqs. (17) or 
(18) (the dotted lines in Fig. 3). On the other hand, in the presence of any randomness defined 
by Eq. (16) in which Ji and J? have the same sign, the quantity | tanh(J/4A;T) | still equals 
unity in the limit T — > 0, and thus this randomness only leads to a quantitative change in the 
thermal dependence of Xxx from that of the uniform case (the dotted lines in Figs, la, Id). In 
Figs. 2 and 3 one can also trace a crossover from a ferromagnetic or antiferromagnetic type 
of thermal behaviour at high temperatures to a paramagnetic type of thermal behaviour as 
T — > 0, which occurs with the introduction of random-bond disorder given by distributions of 
Eqs. (16) with J X J 2 < 0, (17) or (18). Any small amount of disorder causes this crossover, but 
the smaller the disorder the lower is the temperatures at which it occurs. 
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It is also of interest to compare the results for the alternating-bond chains of Section 2 with 
that for the random chain (described by Eq. (16)) having equal concentrations of the two types 
of bond (i.e. c = 0.5). It was seen for the regular alternating-bond models having J x > 
and J 2 < (or vice versa), that the low temperature longitudinal susceptibility either did (in 
the case of J1J1J2J2 ■ ■ ■), or did not (in the case of J\ J2J1J2 ■ ■ •) diverge at low temperatures. 
Further, the divergence when it occurs, is greater than 1/T. For the corresponding random 
chain (c = 0.5) having these signs for the exchange interactions, it has been observed above 
that a 1/T dependence results instead. (Compare the short-dashed and long-dashed lines with 
the curve for c = 0.5 in Figs, lb, lc and 2). On the other hand, when J 1 and J 2 are either both 
positive or negative, the corresponding curves only differ quantitatively (Figs, la and Id). 

Finally, it is noted that the initial longitudinal static susceptibility for chains with Gaus- 
sian and Lorentzian disorders having J = exhibit the temperature dependence of an ideal 
paramagnetic Xxx = 1/4/cT, independent of the strength of the disorder. 

Consider now the thermal behaviour of the initial transverse static susceptibility xJz, Eq. 
(15), (Figs. 4-7). For the uniform Ising chain it is finite at all temperatures; equal to 1/(2 J) 
at T = 0, slightly increasing with increase of temperature, and then decreasing to approach 
the temperature dependence characteristic of an ideal paramagnet at high temperatures. The 
weaker the exchange interaction between the longitudinal (x) components of the spins, the larger 
is the value of Xzz at zero temperature, the smaller is the temperature at which Xzz exhibits its 
maximum, and the smaller is the temperature at which it exhibits behaviour characteristic of an 
ideal paramagnet. These features will be reflected in the sequence of curves for the temperature 
dependence of Xzz of the random chain, described by Eq. (16), as the concentration is varied. 
Fig. 4 shows such a series of curves as c varies from 1 to for the particular case when 
J\ —\ J \ and J 2 = 0.3 | J |. The susceptibilities of some regular alternating-bond chains 
are also depicted. It is interesting to note that the initial transverse static susceptibility for 
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the random-bond chain having c = 0.5 is enhanced over that for the alternating-bond chain 
J\JiJ\Ji--- , coincides with that for the chain J1J1J2J2J1J1J2J2 ■ ■ ■ , and is suppressed in 
comparison to that for the chain J x Jx ^ J2J2J2 ■ ■ ■ , etc.. 

Figs. 5,6 illustrate a difference in the influence of Gaussian, Eq. (17), and Lorentzian, Eq. 
(18), bond disorder on the temperature dependence of xlz- For increasing, but small, Gaussian 
disorder strength an increasing enhancement of Xzl occurs at low temperatures (curves 2-4 
in Fig. 5). Above a certain disorder strength, this low temperature enhancement is reduced 
(curves 5,6), and finally a suppression of \zz below its non-disordered value occurs (curve 7). 
In contrast to this behaviour for Gaussian disorder, chains with Lorentzian disorder exhibit a 
decrease in Xzz with increasing disorder strength for almost all temperatures (Fig. 6). 

Finally, it is noted that in contrast to the initial longitudinal static susceptibility, the initial 
transverse static susceptibility for chains having Gaussian or Lorentzian disorder, with J = 0, 
decreases with increase of disorder strength (Fig. 7). This is not surprising, since the larger 
the exchange coupling, the smaller is Xzz, and with increasing a 2 or F larger values of exchange 
couplings are present more often. 

4 Discussion 

To summarise. An alternative, but more natural derivation than that of Ref. ||, of the exact 
formulae for the initial static susceptibility of random-bond Ising chains has been presented. 
The former requires evaluation of the thermal averaged magnetization to first-order in the 
magnetic field strength, whereas the latter, requires evaluation of just the field-independent 
part of the double-time pair correlation function. The latter approach also permits one to 
derive similarly the frequency dependent susceptibility and the structure factor of nonuniform 
and random spin-| Ising chains; work in this direction is in progress. Dependent on the type of 
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randomness present, it has been demonstrated that the disorder may lead to either a qualitative, 
or to only a quantitative change, in the thermal dependence of the initial static susceptibilities. 
Even a small amount of disorder may lead to a qualitative deviation from that of the nonrandom 
case at low temperatures. A difference in the temperature dependence of the initial static 
susceptibilities of regular alternating-bond chains having two types of exchange bonds, and the 
corresponding random-bond chain, has been revealed. Comparison of the numerical results 
for chains having Gaussian bond-disorder and for those having Lorentzian bond-disorder show 
a marked difference in the dependence of the transverse susceptibility on the strength of the 
disorder. Although, the theoretical results observed in our work should prove valuable in 
understanding the effects of disorder on the observable properties of Ising chain materials, there 
are to our knowledge, no experimental results yet available that enable a direct comparison 
between theory and experiment to be made. However, with for example, the synthesis of 
magnetic chain molecular materials becoming a reality, this lack of experimental data may be 
resolved in the future. 

A part of this work was done during the Summer College in Condensed Matter on "Statistical 
Physics of Frustrated Systems" (July - August 1997) at the ICTP, Trieste. One of the authors 
(O.D.) is grateful to the ICTP for support and hospitality. 
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List of figure captions 



Fig. 1. Thermal dependence of the initial longitudinal static susceptibility of disordered- 
bond (described by Eq. (16)) and regular alternating-bond Ising chains when J\ and J 2 are 
respectively, (a) — | J |, —0.3 | J |, (b) — | J |, 0.3 | J |, (c) | J |, —0.3 | J | and (d) 
| J |, 0.3 | J |. Curve 1: uniform chain, c = 1. Curves 2 and 3: regular alternating chains 
J\ J 2 J\ J 2 • • • , J1J1J2J2 ■ ■ ■ , respectively. Curve 4: disordered chain, c = 0.5. Curve 5: uniform 
chain, c = 0. Curve 6: ideal paramagnet, J\ = J 2 = 0. 

Fig. 2. Thermal dependence of kT\^ for disordered-bond (described by Eq. (16)) and regular 
alternating-bond Ising chains when J\ = — \ J \ and J 2 = 0.3 | J |. Curve 1: uniform chain, 
c = 1. Curves 2-8: disordered chains, c = 0.9, 0.8, 0.7, 0.5, 0.3, 0.2 and 0.1, respectively. Curve 
9: uniform chain c = 0. Curves 10 and 11: regular alternating-bond chains J\ J 2 Ji J 2 • • • and 
J\ J\ J 2 J 2 . . . , respectively. Curve 12: ideal paramagnet, J\ = J 2 = 0. 

Fig. 3. Thermal dependence of kT\^ for Ising chains with (a) Gaussian and (b) Lorentzian 
bond disorder. Curves 1 to 3 are for J = - I J I with (a/ \ J |) = or Ff | J |= 0, 
[aj I J |) 2 = 0.5 or T/ I J |= 0.5 and (a / \ J |) 2 = 1 or F / \ J \= 1, as appropriate to (a) and 
(b), respectively. Curves 4 to 6 are for J =| J |, with the same pairs of {a I I J \Y and F/ \J\ 
as for curves 1 to 3, respectively. Curve 7 is for the ideal paramagnet, Jo = 0. 

Fig. 4. Thermal dependence of the initial transverse static susceptibility of disordered-bond 
(described by Eq. (16)) and regular alternating-bond Ising chains when J\ —\ J \ and J 2 = 
0.3 I J |. Curve 1: uniform chain, c = 1. Curves 2 to 4: disordered chains, c = 0.7, 0.5, 0.3, 
respectively. Curve 5: uniform chain c = 0. Curves 6 to 8: regular alternating-bond chains 
Ji J 2 J\ J 2 ■ ■ ■ , Ji Ji J 2 J 2 ■ ■ ■ and Ji Ji Ji J 2 J 2 J 2 • • • , respectively. 

13 



Fig. 5. Thermal dependence of the initial transverse static susceptibility of Ising chains with 
Gaussian bond-disorder, when J =\ J \. Curve 1: no disorder, (a/ \ J |) 2 = 0. Curves 2 to 7: 
(a/ | J |) 2 = 0.1, 0.2, 0.5, 1, 2 and 4, respectively. 

Fig. 6. Thermal dependence of the initial transverse susceptibility of Ising chains with 
Lorentzian bond-disorder, when J =| J |. Curve 1: no disorder, Y/ \ J |= 0. Curves 2-6, 
r/ | J |= 0.02, 0.1, 0.2, 0.5 and 1, respectively. 

Fig. 7. Thermal dependence of the initial transverse susceptibility of Ising chains with (a) Gaus- 
sian and (b) Lorentzian bond-disorder, when J = 0. Curve 1: ideal paramagnet, (a/ \ J |) 2 = 
or Y / | J |= 0, as appropriate in (a) and (b). Curves 2 and 3: (a/ \ J |) 2 = 0.5 or Y / \ J \— 0.5 
and (a/ \ J |) 2 = 1 or Y / \ J |= 1, respectively. 
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